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ABSTRACT 

We have observed an ./V-body/Smoothed Particle Hydrodynamics simulation of a 
Milky Way like barred spiral galaxy. We present a simple method that samples N- 
body model particles into mock Gaia stellar observations and takes into account stellar 
populations, dust extinction and Gaia's science performance estimates. We examine 
the kinematics around a nearby spiral arm at a similar position to the Perseus arm at 
three lines of sight in the disc plane; (l, b ) = (90, 0), (120, 0) and (150,0) degrees. We 
find that the structu re of t he pe culiar kinematics around the co-rotating spiral arm, 
which is found in [Kawata et al.l (j2014bt) , is still visible in the observational data ex¬ 
pected to be produced by Gaia despite the dust extinction and expected observational 
errors of Gaia. These observable kinematic signatures will enable testing whether the 
Perseus arm of the Milky Way is similar to the co-rotating spiral arms commonly seen 
in IV-body simulations. 

Key words: methods: A-body simulations — methods: numerical — galaxies: struc¬ 
ture — galaxies: kinematics and dynamics — The Galaxy: structure 


1 INTRODUCTION 

The spiral features visible in many galaxies have long 
been the subject of debate. Although it has been almost 
a century since t he re solution of the “great debate” of 
IShaplev fc Curtin dl92ll l. when it was argued over whether 
these beautiful spiral structures were nebulae within our 
galaxy or galaxies in their own right, the mechanisms which 
generate them are still uncertain. One of the problems with 
developing a comprehensive theory of spiral arms is the so 
called “winding dilemma". It is known from observations of 
disc galaxies that the stars in the inner region have a higher 
angular velocity than those in the outer region. Therefore 
the spiral structure should “wind up" relatively quickly if 
the spiral arms rotate at the mean rotation velocity of the 
stars (e.g. IWilczvnskil Il896l l. contrary to observations of 
many “grand design” spiral galaxies. A proposed solution 
to the winding dilem ma is given by spiral density wave the¬ 
ory (Lin & Shul ll964f l which treats the spiral structure as a 
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density wave which can rotate rigidly as a feature with a 
constant pattern speed and thus be long lived. 


However, no IV-body simulations have yet been able 
to reproduce these long lived stable spiral arms, de¬ 
spite the increase in computational power an d resolu¬ 
tion which has occur r ed in recent years (e.g. ISellwoodl 
l201ll : iDobbs &; Babal l2014l l . Recent work has shown 
spiral modes and waves which s urvive over multi¬ 
ple rotations d Quille n et al.l l201ll : iRoskar et al.l 120131 : 
ISellwood fe Carlberd 2014 1 while the spiral arm features 
in the stella r mass are sh ort-li ved but recurren t (e.g . 
Sellwood fe Carlberd l~984l: Carlberg fe Freed man! 1985; 


Bottemal 2003: Fuiii et al.1 201 if Grand et al. 2012 


2013; 


Baba et al .1 120131 : iRoca-F&breea et al.ll2013t ID’ Onghia et al.l 
2013) including in galaxies with a central bar (e.g. 
Grand et al.ll2012bll . implying that the large spiral arms vis¬ 
ible in external galaxies may only appear to be rigid struc¬ 
tures extending over the disc, while in fact being made of 
transient reforming features. 


The interpretation of the transient and recurrent spiral 
arm features obse rved in iV- body simu lations is still in de¬ 
bate. For example. iMinchev et al l (120121 ') showed for the first 
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time (by studying the time evolution of the disc power spec¬ 
trum) that spiral wave modes in Wbody simulations can 
last for as long as 1 Gyr, which can justify treating the wave 
modes as quasi-stationary structure, and the transient and 
recurrent spiral arm features can be explained by the super- 
positi on of different modes with different pattern sp eeds (see 
also iRoskar et_aL 2012^ Sellwood&Carlberg[ 2014). On the 
other hand. [Grand et alT ~i 2012^ )7 D’Onghia et ahl d2013h : 
iBaba et al.l 11201311 demonstrated non-linear growth of the 
spiral arm features due to similar but different (i n terms 
of ev olution) mechanisms from swing-amplification llToomrel 
Il98ll) . which could be difficult to explain with the linear su¬ 
perposition of the wave modes. 

Our position within the Milky Way gives us a unique 
view of these spiral structures seen in external galaxies, 
but it comes with its own set of problems which we must 
overcome when studying them. The location and kinemat¬ 
ics of the gaseous component of the ar ms ma y be deter ¬ 
mined from HI and CO observations (e.g. Dame et al.ll200l] : 
iNakanishi fc Sofu3120031 : iKalberla fe Kerdl2009l) . However 
to observe the kinematics of the stellar component in and 
around the spiral arms we must look through the disc plane, 
which carries the heaviest levels of dust and gas, and thus 
high levels of extinction. 

Dust extinction has long been a problem for Milky Way 
model construction. Although there are reasonably reliable 
extinction maps for extra galactic sources whose extinc¬ 
tion by the interstellar medium of t he Milky W ay can b e 
corrected as a function A\(l,b) (e.s. ISchlegel et ahlllOQsl ). 
three dimensional extinction mapping for sources within the 
Milky Way i.e. A\(l,b,d) is more challenging. There are 


three dimensional extinction maps for individual sections of 


the sky (e.g. Drimmel & Spergel 2001; Marshall et al. 

2006; 

Hanson & Bailer-.Tones 201-1; Sale & Magorrian 2014 

) and 


two dimensional map s have been extended to three dimen¬ 
sions (e.g. iDrimmel et al.l [20031) . Howe ver a truly Galacti c 
3D extinction map does not yet exist llRix fc Bovvl f2013f ). 
The European Space Agency (ESA)’s Gaia mission will 
help us map the stellar structure and kinematics of the 
Milky Way, and help con strain extinction at the same time 
llBailer- Jones et al .11 201 .‘3) . 

Gaia, which was launched on the 19th December 
2013 will provide detailed astrometric (e.g. iLindegren et al.l 
l2012ll . s pectroscopic (e.g. iKatz et al.l [201 il l and photomet¬ 
ric (e.g. IJordi et al.l |201y) information for around one bil¬ 
lion stars in the Milky Way. Detailed information o n Gaia 
scient ific accuracies is available in, for example, Ide Bruii~n3 
(I2012H . Synthetic Gaia mock data have already been used 
to demonstrate d ifferent app li cation s of the real Gaia data 
set. For example. lAbedi et al.1 d2014ll use three tracer popu¬ 
lations (OB, A and Red Clump stars) with the Gaia se¬ 
lection function, errors and dust extinction, and demon¬ 
strated that the Gaia mock d ata can recover the p arameters 
of the Galactic warp. iRomero-Gomez et al.l (l2014ll examine 
the Galactic bar in the Gaia observable space using Red 
Clump tracers with the Gaia selection function, errors and 
dust extinction combined with selected Red Clump stars 
from the Apache Point Observat ory Galactic Evo lution Ex¬ 
periment (APOGEE DR10, e.g. lAhn et ahl l2014h showing 
th e value of combining d ata from complimentary surveys. 
In lHunt fe Kawatal (120141 ) we show that we can recover the 
large scale structure of the Galactic disc with our Made-to- 


Measure Galaxy modelling c ode, primal Jllunt et al.ll20l.j ; 
iHunt fe Kawatal I2013I . I2014I) . and make a good estimation 
of the patten speed of the bar, using tracer populations of 
M0III and Red Clump stars with the Gaia selection func¬ 
tion, errors and dust extinction. 

There exist full mock catalogues of Gaia stars, e.g. the 
Gaia Universe Model Snapshot (gums) which provides a 
view of the Besan qon Galaxy model as seen from Gaia 
l|Robin et al.l 120121 ) . taking into account dust extinction 
while assuming there are no observational errors. This de¬ 
tailed prediction of Gaia observations gives an excellent in¬ 
dication of the volume and quality of data which will be¬ 
come available from Gaia , predicting 1.1 billion observable 
stars, almost 10,000 times more than from its predecessor 
Hipparcos. gums can be extended th rough the Gaia Ob¬ 
ject Generator (gog) liLuri et al . 20H) to simulate inter¬ 
mediate and final catalogue data including the introduction 
of realistic astrometric, photometric and spectroscopic ob¬ 
servational errors to the catalogue based upon Gaia science 
performance estimates. While these mock data provide an 
excellent example of the capabilities of Gaia, the Besanqon 
galaxy model is an axisymmetric model and a kinematic 
model not a dynamical model. Although Gaia will not pro¬ 
vide accelerations, the kinematics it will provide are from 
a dynamical system, the Milky Way. Thus it is important 
for our purpose to generate catalogues from fully dynami¬ 
cal models with non-axisymmetric structures, such as spiral 
arms and a bar, which for example IV-body disc galaxy mod¬ 
els can provide. 

Therefore we propose here to create mock Gaia obser¬ 
vations from an ./V-body m odel us ing a population synthesis 
code such as GALAXlAjShmrm^etffi. 201l|), o r the method¬ 
ology presented in IPasetto et al.l 1 20121) or iLowing et al.l 
|20 i 4 GALAXIA is a flexible population synthesis code for 
generating a synthetic stellar catalogue from an ./V-body 
or an analytical galaxy model over wide sections of the 
sky, with a sampling scheme which generates a smoothly 
distributed sample of stars. Synthetic catalogues generated 
from dynamical Galaxy models are essential for preparing 
to exploit the real Gaia catalogue and can be used to deter¬ 
mine whether certain features within the Milky Way will be 
visible to Gaia. 

In our previous work dKawata et al.l [2014bl) we exam¬ 
ined the kinematics of both the stellar and gas components 
around a transient, co-rotating spiral arm in a simulated 
barred spiral galaxy similar in size to the Milky Way. Al¬ 
though this arm is transient, similar arms recur during the 
evolution of the galaxy. We made predictions of observable 
kinematic signatures that may be visible in the Milky Way’s 
Perseus arm if it is also a transient, recurrent and co-rotating 
spiral arm. We then compared our simulation with dat a from 
APOGEE and the maser sources from iReid et ahl (|2014l ) 
measured by the Bar and Spiral Structure Legacy (BeSSeL) 
survey and the Japanese VLBI Exploration of Radio As¬ 
tronomy (VERA), finding tentative agreement between our 
simulation and the observations. Owing to the low number 
of maser sources and the lack of distance information for the 
APOGEE stars no firm conclusions could be drawn; however 
it is encouraging to see similar features in both, including 
the possible signatures of a co-rotating spiral arm. 

In this paper we build upon the previous work by gen¬ 
erating a stellar sample with different populations from the 
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simulation data in iKawata et al.l (l2014bh and making mock 
observations of these stars taking into account the expected 
Gaia science performance estimates. The aim is not to make 
further predictions about the kinematics of transient, re¬ 
current and co-rotating spiral arms but rather to examine 
whether these signatures, remain visible in the Gaia data if 
they exist in the Milky Way. 


2 SIMULATION 


We use the simula ted galaxy which is p resented in 
IKawata et al.l (l2014bh and lGrand et al.l d2014bf >. The details 
of the nume rical simulation code, and the galaxy model are 
described in iKawata et al.1 J2014o ). We briefly describe the 
galaxy model in this section. The galaxy is set up in isolated 
conditions, and consists of a gas and stellar disc but no bulge 
component. The di scs are embedded in a static dark mat¬ 
ter ha lo potential llRahimi fe Kawatal 120121 : IKawata et al.l 
l2014bl) . The dark matter halo mass is Md m = 2.5 x 10 12 
M@, an d the dark matter de nsity follows the density pro¬ 
file from lNavarro et al.l (|l997f ). with a concentration param¬ 
eter of c = 10. The stellar disc is assumed to follow an 
exponential surface density profile with the initial mass of 
Md,* = 4.0 x 10 10 Mg, a radial scale length of lid,* = 2.5 
kpc and a scale height of zg ,* = 350 pc . Th e gas disc is set 
up following the method of ISpringel et al.1 ( 20050 . and has 
an exponential surface density profile with the scale length 
of Rd, g = 8.0 kpc. The total mass of the gas is 10 10 Mg. 
The simulation comprises 10 6 gas particles and 4 x 10 6 star 
particles; therefore each particle has a mass of 10 4 Mg. The 
resolution is s ufficient to minimise numerical heating from 
Poisson noise dFuiii et al.ll201ll : ISellwoodll2013l ). We employ 
a minimum softening length of 158 pc (equivalent to a Plum¬ 
mer softening length of 53 pc) with the spline softening and 
variable softening length for gas particles as suggested by 
IPrice fe Monaghanl 120071 ). 

The radial profile of the mean metallicity of stars and 
gas is initially set by [Fe/H](I?) = 0.2 — 0.05(12/1 kpc), and 
the metallicity distribution function at each radius is centred 
on the mean metallicity value with the dispersion set to a 
Gaussian distribution of 0.05 dex for the gas and 0.2 dex for 
the stars. The stellar ages are set randomly between 0 and 
10 Gyr for stars present at the beginning of the simulation. 

The simulation was run for 1 Gyr from the ini¬ 
tial conditions with the IV-bod y smoothed partic l e hy¬ 
drodynamics code, gcd+ (e.g. Kawat^&Gibson 2003L 
Rahimi fc Kawatal 120121 : iBarnes et al.l 120121 : IKawata et al.l 


20131 . 2014al ) without the inclusion of any continuous ex¬ 


ternal inflow of gas for simplicity. In this paper we use the 
same snapshot of the galaxy as used in lKawata et al.l (|2014bi ) 
which is taken at t = 0.925 Gyr, as this snapshot shows a 
spiral arm at a similar location to the Perseus arm from the 
Milky Way (see Fig. [TJ. 


3 GAIA MOCK CATALOGUE 

In IKawata et al.l (l2014bh the kinematics of the spiral arm 
shown in Fig. [T]are examined at three lines of sight Ilos = 
90, 120 and 150 deg, with &i os = 0 because of the lower ex¬ 
tinction relative to other lines of sight in the plane. Predic¬ 



Figure 1. Snapshot of the simulated galaxy from K awata et al. 
I2014bl ) which is also used in this paper. The left (right) panel 
shows the face-on view of the star (gas) particle distribution. The 
solid line indicates the position of the spiral arm identified. The 
observer is assumed to be located at (x,y) = (—8,0) kpc. Three 
line-of-sight directions (Ilos = 90,120 and 150 deg) are high¬ 
lighted with the dotted lines. The galaxy is rotating clockwise. 


tions are made of the observational signatures of co-rotating 
spiral arms notably the difference in kinematic structure be¬ 
tween the trailing near side and leading fa r side of the spiral 
arm. In gene r al. inlKawata et all d2014bl ) (as also shown in 
iGrand et al.l ll2014ah ) the stars in the trailing near side ro¬ 
tate slower because they tend to be at the apo-centre and 
migrate outward, and the stars in the leading far side ro¬ 
tate faster as they tend to be at the peri-centre and migrate 
inward. There are however some stars which follow the op¬ 
posite trend, leading to multiple populations seen in the ro¬ 
tational velocity in the leading far side; one faster, and one 
slower than the single population in the trailing near side. 
These features which will be discussed later may be caused 
by the co-rotation resonance of the spiral arm, and are vis¬ 
ible at different galactic longitudes because the arm in the 
simulati on co-rotates at all the examined radial range. How¬ 
ever, in lKawata et al.l (l2014bl ). the spiral arm kinematics are 
examined using the full, error and extinction free IV-body 
data and thus such trends, when present, are easy to iden¬ 
tify. 

In this Section we describe how we generat e a sam ple 
of stars from the IV-body model of lKawata et al.l (l2014bl) to 
produce a mock Gaia catalogue. It is worth noting th at the 
population synthesis code, galaxia aSharma et al.l 1201 In 
provides a tool to generate stellar populations from IV-body 
simulation data. However, because we plan to combine such 
a tool with our Made-to-Measure Galaxy modelling code, 
primal, we have developed our own simplified version of 
galaxia, a population synthesis code called snapdragons, 
(Stellar Numbers And Parameters Determined Routinely 
And Generated Observing IV-body Systems), snapdragons 
uses the same isochrones and extinction map as galaxia, 
but uses a different and more simplistic process to gener¬ 
ate the stellar catalogue which is described in Section 13.21 
snapdragons allows us to add the expected Gaia errors 
more easily, and enables us to track the link between sam¬ 
pled stars and their parent IV-body particle for our future 
studies, e.g. primal modelling of the Galactic disc by fitting 
tracers from multiple stellar populations, and identifying ra¬ 
dially migrating stars and non-mig rating stars trapped by 
the spiral arm llGrand et al.ll2014al ). 
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3.1 Extinction 


We use th e extinction map o f the Milky Way taken from 
galaxia (|Sharma et al.11201 11 1, which is a 3D polar loga¬ 
rithmic grid of the du st extinction constructed u sing the 
method presente d in Bland-Haw thorn et al. I I 20 T 0 I) and the 
dust maps from ISchleeel et al.l (1 1998T) . The same extinc¬ 
tion is applied in iHunt fc Kawatal (|2014 ) and more detai l 
is given there. In an update from IHunt fe Kawatal (l2014l l 
we fo llow the correc tion to the Schlegel Eb-v presented in 
ISharma et al. Mioll) such that 


Eb-v = Eb-v ^ 0.6 + 0 . 2^1 — tanh ^ ^ 8- ^ — °' 15 

(1) 

This correction is made as it has been suggested (e.g. 
lArce &: Goodman] 1 19991 : lYasuda et al.l l2007li that the red¬ 
dening is overestimated by the maps from ISchlegel et al.1 
(Il998h by ~1.3-1.5 in regions with high extinction with 
Av > 0.5 ( Eb-v > 0.15). This correction reduces extinc¬ 
tion by ~ 40% for low latitude high extinction regions but 
has minimal effect on high latitude low extinction regions. 



3.2 Population Synthesis: snapdragons 


The goal of this population synthesis code is to split each N- 
body particle from the galaxy simulation into an appropriate 
number of stellar particles creating a mock catalogue of ob¬ 
servable stars from our A-body model. We must choose an 
IMF and a set of isochrone s with which to work. We choose 
a Salpeter IMF (|Salpeteij 1 19551 ) where the IMF, 4>(m), is 
defined in each mass interval dm as 

$(m)dm = Am~^ x+1 ^ Am, (2) 


where x = 1.35 is the Salpeter index, and A is a constant 
for normalisation in the desired mass range. We set this 
constant as 



— X 


m 



-1 


(3) 


where mi is the A-body particle mass, m*,;, max is the max¬ 
imum initial mass of any surviving star and m* !m i n is the 
minimum stellar mass to be considered. We make use of 
the P adova isochrones fe.g. lBertelli et al . 1994; Marig o~et al.1 
2008), although the choice of isochrones (and IMF) may be 
substituted with others with no change to the methodology. 

It is worth noting that the Padova isochrones are avail¬ 
able only for stellar masses abov e 0.15 Afp. ga l axia for 
example uses the isochrones from IChabrier et al.l (|2000T) to 
extend the mass limit down to 0.07 Mg, which is the hy¬ 
drogen mass burning limit. We set our lower limit on stellar 
mass as m«, m j n = 0.1 Mg to c orrespond with the simula¬ 
tion from [Kawata et al.l ( 2014bf ) and extrapolate from the 
Padova isochrones for 0.1 < Mq < 0.15. It is relatively safe 
to do this because all such stars lie on the main sequence. 
Additionally these exceedingly faint stars will not be visible 
at the distance of the spiral arms which are the focus of this 
work. 

As discussed in Section [5] each A-body star particle in 
the simulated galaxy has been assigned an age and metallic- 
ity within the chemodynamical code GCD-f , then it is made 
to evolve. When we examine the snapshot, each particle is 
matched to its nearest isochrone in both metallicity and 


age from the grid of isochrones which are extracted from 
galaxia. Once an isochrone is selected, we identify m*,j, max 
from the isochrone. We then determine how many stars to 
sample from the A-body particle by integrating the IMF 
over the desired mass range; 

max . 

N s = A m~ 1+1 dm, (4) 

where m* il! <i/ lirn is minimum mass required for the star par¬ 
ticle to be brighter than our apparent magnitude selection 
limit, Viim, taking into account the extinction value at the 
position of the parent particle. Stars smaller than 
are not used in the subsequent analysis, to save on compu¬ 
tational time. 

We then randomly sample stellar masses from the sec¬ 
tion of the isochrone N s times. We have weighted the ran¬ 
dom selection by the IMF using the equation 

m* = (Rm“* max + (1 - R)m~l <Vu J ^, (5) 

where R is a random number between 0 and 1. 

The isochrones are comprised of discrete stellar data, 
and therefore we then interpolate within the nearest 
isochrone values of Mv and V — I c to determine My* and 
V — I c * for the generated m*. At this stage we assume the 
generated stars have the same position and velocity as their 
parent particles. 


3.3 Observational Errors 

Having generated the visible stellar catalogue we then add 
observational errors based upon the Gaia Science Perfor¬ 
mance estimate^]. We use the post launch error estimates 
approximated from the estimat es in pre-launch perfor mance 
by Merce Romero-Gomez fe.g. IRomero-Gomez et al.ll2014il . 
provided through the Gaia Challenge collaboratiorQ. We as¬ 
sume the position and velocity of the Sun is known. We 
locate the observer at (—8,0,0) kpc as shown in Fig. [l] and 
the motion of the Sun is assumed to be 228 km s -1 . For this 
work, while generating the stellar catalogue we produced 
stars only brighter than V|j m < 16 mag, which is well within 
Gaia's rriG < 20 mag magnitude limit for the astrometry. 
However, because we are interested in the Galactic radial 
and rotation velocity for the stars, which requires the full 
6D phase space information, we chose the lower magnitude 
limit where Gaia RVS can produce the reasonably accurate 
line-of-sight velocity. Note that the errors are added to the 
parallax, proper motion and line-of-sight velocities. 

A full description of the method to add the pre-launch 
Gaia error is available in iHunt fc Kawatal (120141 4. However 
the Gaia science performance estimates have been revised 
after launch, and as such a correction must be made. The 
error in parallax has increased, and although it has little 
effect for stars with mv < 16 mag which we work with in this 
paper, the coefficients within the equation to describe the 
pre-launch parallax performance (provided by Kazi, Antoja 
& DeBruijne (Oct. 2014) by fitting to the new estimations 


1 http://www.cosmos.esa.int/web/Gaia/science-performance 

2 http://astrowiki.ph.surrey.ac.uk/dokuwiki/doku.php 
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on the Gaia science performance web page) are revised to 

(j n = (—11.5 + 706.1z + 32.6 z 2 ) 1/2 

x (0.986 + (1 - 0.986)(U - I c )), (6) 

where 

z — max(10°' 4(12_15 \ io 0 ' 4< - G-15) ), (7) 

correcting also the t ypo for equations © and © in 
[Hunt fc Kawatal 120141 ). 

Additionally, because of the loss of spectroscopic accu¬ 
racy by ~ 1.5 mag in the RVS post launch performance we 
also apply a correction to the error function for the end of 
mission radial velocity. We change the tablcjf] of values a and 
b, again determined by fitting the revised performance es¬ 
timates on the Gaia science performance web page, for the 
equation 

a Vr = l + be a(v ~ 14) , (8) 

where a and b are constants dependant on the spectral type 
of the star. The new table along with the code to add the 
Gaia error is available onlincQ 


4 RESULTS 

As disc ussed in Section [3] it was shown in iKawata et al.l 
(l2014bl l that in general the stars in the trailing near side of 
the spiral arm rotate slower than average because they tend 
to be at the apo-centre, and the stars in the leading far side 
of the spiral arm rotate faster than average as they tend 
to be at the peri-centre. However, there are groups of stars 
which follow different trends leading to multiple populations 
which will be discussed later. It is important to determine 
whether such features will still be visible in the Gaia cata¬ 
logue, not just the error and extinction-free A-body model. 
In this Section we show the result of sampling these A-body 
data into stellar data, first looking at the properties of the 
resulting mock stellar catalogue, and then examining the spi¬ 
ral arm kinematics with the stellar data taking into account 
dust extinction and Gaia science performance estimates. 

4.1 Population synthesis 

In this section we describe the stellar catalogue produced 
by snapdragons, and show the resulting colour magnitude 
diagram (CMD) varying the area of the sky coverage. Fig. [2] 
shows the CMD for stars generated by snapdragons from 
particles within a square region of ±2 deg (upper) and ±5 
deg (lower) around ( l , b) = (90, 0) deg. The upper panel 
of Fig. [2] shows clearly the individual stellar isochrones be¬ 
cause there are only a small number of A-body particles in 
the selected region, and each particle has only one age and 
metallicity. These problems are resolved when smoothing is 
applied in the p hase space di s tribut ion and age-metallicity 
distribution fe.g. lSharma et ahll201lll . However, as discussed 
in Section ^.2l we deliberately avoid this smoothing to main¬ 
tain the clear particle-star relation. The lower panel of Fig. 
[2] shows no such discrete structure, as there are sufficiently 

3 http://www.cosmos.esa.int/web/Gaia/table-5 

4 https://github.com/mromerog/Gaia-errors 



V-I c 


Figure 2. Colour magnitude diagram for stars generated by 
snapdragons from particles within a square region of ±2 deg 
(upper) and ±5 deg (lower) around ( l,b ) = (90,0) deg. Stars 
with apparent magnitude of my < 16 only are included. 

many particles to cover a broad range of stellar ages and 
metallicities in the CMD. Therefore, care is required with 
the resolution of the A-body simulation and the selection 
function if we discuss in detail the stellar population distri¬ 
bution in the CMD. However, this is unlikely to affect the 
study in this paper. 

4.2 Observable Spiral Arm Kinematics 

In this section we examine if the possible kinematic signa¬ 
tures of co-rotating transient an d recurrent spiral arms iden¬ 
tified in iKawata et al.l (l2014bl) will be visible in the Gaia 
data even given the dust extinction in the disc and Gaia's 
science performance accuracy. A d etailed analysis of th e 
kinematics themselves is the focus of lKawata et al.l (l2014bl) . 
while this work is concerned with the visibility of this kine¬ 
matic structure in the Gaia data. We examine the rotational 
velocities o f the stars in the catal ogue for different distances 
because in iKawata et al. . ( 2014b ) we found the rotation ve¬ 
locity is most affected by the transient co-rotating spiral 
arm. Then we calculated the Probability Density Function 
(PDF) of the rotation velocity of stars behind and in front 
of the spiral arm using Kernel Density Estimation (KDE) 
which we are using as a desirable alternative to histograms 
(e.g. 1 Wassermanl 2()() <i i. 

Fig. EH shows a smoothed contour plot of the galacto- 
centric rotational velocity against distance for particles and 
stars within a square region of ±5 degrees around (l, b) = 
(90,0) (left), (l,b) = (120,0) (middle) and ( l,b ) = (150,0) 
(right). This compares the kinematics of the underlying A- 
body model (upper) with the stellar catalogue generated 
with snapdragons, before (middle) and after (lower) the 
addition of the errors from the Gaia science performance 
estimates. Owing to the high percentage of low mass and lu¬ 
minosity stellar types which would dominate the selected re¬ 
gion and saturate the plot at small distances, we have made 
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cuts to our sample to visualise the underlying kinematic 
structure from the stellar catalogue. We have first cut the 
sample of stars in all three lines of sight with absolute mag¬ 
nitude, My < — 1, calculated from the apparent magnitude 
my and observed distance d 0 b s , assuming the dust extinc¬ 
tion is known. We then cut with cr„ los / (vi os x dobs) < 0.015 
kpc -1 to select the stars with lower error in the line of sight 
velocities at a smaller distance to generate similar quantities 
of data at different distance scales. This is purely for illus¬ 
tration purposes and we are not suggesting that this is the 
best possible selection function. The upper panels of Fig. [3] 
show the different kinematic structure in the IV-body model 
at the different lines of sight. These are the same data as 
t hose s hown in the top panels of Fig. 4 from iKawata et all 
( 2014bh . Note that the density colour scale for the ZV-body 
data is different from the stellar data in the middle and lower 
panels. 

The middle row of panels of Fig.[3]show the velocities of 
the selected stars, which appear slightly different from those 
of the IV-body data owing to the selection function. While 
the general shape of the distribution has been recovered, 
at (Z, b ) = (90, 0) deg (middle left) the fast rotating stars 
within the arm dominate the density scale and wash out the 
rest of the plot slightly. At (Z,6) = (120,0) deg (middle), 
although there is some saturation around 220 km s -1 the 
kinematic structure is clearly visible and is a good match to 
the particle data. Similarly at ( l,b ) = (150,0) deg (middle 
right), despite the lower number of counts, the kinematic 
structure is clearly shown. 

The lower panels of Fig. [3]show the error affected rota¬ 
tion velocity and distance for the selected stars taking Gaia 
science performance estimates into account. The rotation ve¬ 
locity is calculated from the observed parallax, proper mo¬ 
tion and line of sight velocities. At (l, b) = (90,0) (lower 
left) the shape of the distribution remains relatively un¬ 
changed, with the main loss in accuracy occurring around 
dobs ~ 7 — 10 kpc. The recovery of the kinematic struc¬ 
ture around the spiral arm around d 0 bs ~ 4 kpc remains 
almost identical to the case without observational errors. At 
(l,b) = (120,0) (lower middle) the visible loss of accuracy 
is again in the outer region of d 0 bs ~ 7 — 10 kpc, with the 
region containing the spiral arm remaining very similar to 
that of the error free case. At (l,b) = (150,0) (lower right), 
the entire distribution remains very similar to the middle 
right panel, the case without Gaia like observational errors. 

Fig-EJshows the PDF’s, with a KDE bandwidth of 4, for 
the rotational velocity of the stars in the catalogue within 
a square region of ±5 degrees around (l,b) = (90,0) (left), 
(Z,fe) = (120,0) (middle) and {l,b) = (150,0) (right) in the 
trailing near side, between 1 and 2 kpc closer than the cen¬ 
tre of the arm (upper) and leading far side between 1 and 
2 kpc further than the centre of the arm (lower). Note that 
these distance bins were chosen as they show the discussed 
structure most clearly; the same features are present closer 
to the arm but are less clear. The centre of the arm was de¬ 
termined to be at d = 4.0 kpc at (Z, b) = (90, 0), d = 3.4 kpc 
at (Z, b) = (120, 0) and d = 3.3 kpc at (Z, b) = (150, 0). Note 
that Fig. [4] uses all the stars with my < 16, not applying 
the selection function used for illustration purposes in Fig. 
[3] At all three lines of sight Fig. [4]shows a clear difference in 
the distribution of velocities for the ‘true’ data (black solid) 
when comparing the different observed distances, as shown 


in lKawata et al.l l2014bh . This is a positive outcome consid¬ 
ering the loss of data from the dust extinction. When com¬ 
paring the ‘true’ (black solid) stellar catalogue data with the 
stellar data taking into account dust extinction and Gaia's 
expected errors (red dashed) a general smoothing out of the 
structure is evident in the ‘observed’ data. The upper panels 
of Fig. [4] showing the trailing near side of the arm show very 
similar PDF’s when comparing the true and observed stel¬ 
lar data, whereas the lower panels showing the leading far 
side show an information loss, especially at (Z,fe) = (90,0), 
where the three peaks are no longer resolved. This is to be 
expected because of the higher distances and therefore ad¬ 
ditional extinction; however at ( l,b ) = (120,0) and (150,0) 
even on the far side of the spiral arm the structure within 
the distribution is still clearly visible. 

When comparing the ‘observed’ data in Fig. [4] in front 
and behind the spiral arms, we see a clear difference in the 
PDF at all three lines of sight. In each case, the PDF in 
the trailing near side of the spiral arm forms a single cen¬ 
tral peak similar to the mean rotation velocity, with a small 
tail towards faster rotation velocities whereas the leading 
far side of the spiral arm shows a broader distribution of 
velocities with a peak velocity faster than the peak for the 
trailing near side. The difference is particularly apparent at 
(Z,fe) = (120,0) deg where the leading far side shows two 
clear peaks, one faster and one slower than the single peak 
in the trailing near side. This bimodal distribution can also 
be seen in the lower middle panel of Fig. [3]between 4.39 and 
5.39 kpc (although note that Fig. |3]uses a different selection 
function). Also at (Z, 6) = (150, 0) deg the single broad peak 
in the trailing near side is easily distinguishable from the 
leading far side which shows three peaks. These three peaks 
are also partially visible in the lower right panel of Fig. [3] 
between 4.2 9 and 5.29 kpc . These features all match those 
observed in iKawata et al. (|2014b ) despite the addition of 
dust extinction and observational errors to the data. 

In general, as shown in lGrand et al.l ll2014al f. the stars 
in the leading side rotate faster as they tend to be at peri- 
centre phase and migrating inward, and stars in the trailing 
side rotate slower as they tend to be at apo-centre phase 
and migrating outward. This explains the single large peak 
in the trailing side, and the largest peak on the leading side 
which has a higher rotational velocity than the single peak 
on the trailing side as shown in Fig. [4] However, when the 
transient spiral arm starts forming, stars which are close to 
the arm on the trailing side and are close to the peri-centre 
phase, are accelerated towards the arm, passing through and 
then slowing down as th ey reach the apo-centr e on the lead¬ 
ing side as discussed in IKawata et al.l (l2014bll . These stars 
correspond to the ‘slower’ peaks visible in the lower panels 
of Fig. [4] Similarly, the stars which are close to the arm and 
close to the apo-centre phase on the leading side are decel¬ 
erated by the arm, and are overtaken by the arm. Then they 
are accelerated again by the arm once they are on the trail¬ 
ing side at peri-centre phase, which corresponds to the small 
tail present at high velocities in the upper panels of Fig. [4] 
The difference in the rotation velocity distribution between 
the leading and trailing side of the spiral arm seen in Figs. 
03 and Q] is that the latter population is smaller than the 
former. It appears that it is easier for stars to escape from 
the arm on the leading side than the trailing side. From our 
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Figure 3. Smoothed linear scale contour plot of galactocentric rotation velocity of simulation particles (upper), selected snapdragons 
stars (middle) and selected snapdragons stars observed with Gaia error (lower) for ( l,b ) = (90,0) (left), ( l,b ) = (120,0) (middle) and 
(/, b) = (150, 0) (right). For the snapdragons stars (middle and lower panels), a limited selection of My < — 1 calculated using my and 
dobs and assuming a known extinction, along with a Vr /(v r x d 0 bs) < 0.15 is shown to avoid overly dense populations of fainter stars at 
smaller distances. This is to visualise the data set, and these faint stars contribute to the subsequent analysis. Note while consistent for 
snapdragons stars across the different lines of sight, the density scale is different for the simulation particles; however the choice of the 
scale is arbitrary. 


analysis of ./V-body simulations this appears to be a common 
feature of transient and co-rotating spiral arms. 

IComparetta &; Quillenl (|2012h propose that the radial 
overlap of multiple longer-lived patterns moving at different 
pattern speeds can reproduce the transient spiral features, 
which when strong enough can lead to radial migration away 
from the co-rotation radius associa ted with co-ro t ating s pi- 
ral arms as seen, for example in iGrand et al.l (|2012al lbh. 
In such a scenario, the spiral arm features are co-rotating, 
which may give rise to the co-existence of many inner and 
outer Lindblad resonances in a range of radii and lead to the 
features visible in Figs. [3] and [4] However, further analysis 
of the spiral arms in IV-body simulations is required before 
drawing firm conclusions on the mechanism that generates 
such kinematic signatures, which we will tackle in future 
studies. 

From Figs. [3] and [4] we find that Gaia's scientific accu¬ 
racy ought to be sufficient to examine the kinematic struc¬ 
ture of the nearby spiral arms in the Milky Way, even on 
the far side of the arm. Fig. [4] shows clear differences in 


the kinematics in the leading and trailing sides of the spiral 
arm, notably the difference in the number and locations of 
the peaks, and the small high velocity tail present in the 
trailing near side. The comparison between the middle and 
lower panels of Fig. 0 shows little difference, implying that 
the observational error from Gaia will have limited effect on 
our ability to study the kinematics of the spiral arms. Fur¬ 
ther examination of galaxy models constructed using the 
different theories of spiral arm formation will be essential to 
determine the distinct kinematic signatures of each theory. 


5 SUMMARY 

We observed our TV-body/SPH simulation of a Milky Way 
like barred spiral galaxy to create a mock Gaia stellar cat¬ 
alogue, with particular interest in the stellar kinematics in 
and around the spiral arms. We foc used on t he sa me thre e 
lines of sight in the disc plane as iKawata et al.l (l2014blf . 
( l,b ) = (90, 0), (120, 0) and (150,0) deg and analysed the 
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Figure 4. Comparison of the distribution of galactocentric rotational velocities for the stars generated by snapdragons within a square 
region of ±5 degrees around (Z, 6) = (90,0) (left), (l,b) = (120,0) (middle) and (1. b) = (150,0) (right) in the trailing near side (upper) 
and leading far side (lower) of the spiral arm which meet the my < 16 selection limit. The black solid curve shows the true velocities, 
and the red dashed curve shows the distribution once the Gaia errors have been applied. The vertical lines show the circular velocity 
(dotted) and the mean rotation velocity (dash-dotted) at the radius of the spiral arm. 


galactocentric rotational and line of sight velocities of the 
selected stars as a function of the distance from the ob¬ 
server. In agreement with existing literature on IV-body spi¬ 
ral galaxy simulations the spiral arm features seen in the 
stellar mass in our model are transient, recurrent and co¬ 
rotating, i.e. the spiral arm is rotating at the circular veloc¬ 
ity of the stars at the selected lines of sight. 


We sh ow that the stru cture in the kinematics identi¬ 
fied in lKawata et al.l ( 2014bi ) remains visible after the inclu¬ 
sion of dust extinction and observational errors based upon 
Gaia science performance estimates. Although the inclusion 
of these observational effects makes the trends less clear, 
they are still observable in the mock Gaia data in front of, 
inside and behind the spiral arm. The structure on the trail¬ 
ing near side is relatively unchanged, whereas the structure 
on the leading far side is, unsurprisingly, more affected, al¬ 
though the bi-modal (or more) and broader distribution of 
the rotation velocities is still clearly visible. Because we be¬ 
lieve that these kinematic signatures are indications of tran¬ 
sient and co-rotating spiral arms owing to the co-rotation 
resonance at all radii, we predict they should be visible in 
the Gaia data at different longitudes if the Milky Way’s 
Perseus arm is also a transient and co-rotating spiral arm. 


Encouraged by the success of this study, we intend 
to repeat the analysis with simulated galaxies which use 
different theories of spiral struct ure formation, for exam¬ 
ple test particle simulations (e.g. Minchgv_fc_Quillgn 2008); 
MmchevetjdJ 201(1 Minchev fc Famaevll201ol : Faur e et alj 


20141 : Antoia et al.l 1201 41 1 and iV -body simulations with a 


fixed spiral arm potential (e.g. IWada et al.l [20111 ). From 
these analyses we expect to make predictions of the kine¬ 
matic signatures of different spiral arm theories, which can 
be tested by the Gaia stellar catalogue. 
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